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We report on-off intermittency in electroconvection of nematic liquid crystals driven by a dichoto- 
mous stochastic electric voltage. With increasing voltage amplitude we observe laminar phases of 
undistorted director state interrupted by shorter bursts of spatially regular stripes. Near a critical 
value of the amplitude the distribution of the duration of laminar phases is governed over several 
decades by a power law with exponent —3/2. The experimental findings agree with simulations of 
the linearized electrohydrodynamic equations near the sample stability threshold. 

PACS numbers: 05.40. +j, 47.20.-k, 47.54. +r, 61.30-v 



a 
'^ 



> 

O 

o 

(N 

O 
0\ 
0\ 

'o 



Oh- 



% 



Systems at a threshold of stability driven by a stochas- 
tic or chaotic process coupling multiplicatively to the sys- 
tem variables may exhibit on-off intermittency character- 
ized by specific statistical properties of the intermittent 
signal. Quiescent (or laminar) periods (off-states) are in- 
terrupted by bursts of large variation (on-states) ; the du- 
ration of laminar periods is governed by power laws with 
exponents universal over a broad class of different sys- 
tems. Early theoretical studies considered systems with 
few degrees of freedom modeled by differential equations 
[yi and mappings 0. There is increasing interest in sys- 
tems with many degrees of freedom ||^ , described by ran- 
dom map lattices Q , larger systems of coupled nonlinear 
elements [g|, and partial differential equations |g,||. Ex- 
perimental results are available mainly for nonlinear elec- 
tric circuits |7|] ; on-off intermittency was also observed in 
a spin wave experiment [p|, in optical feedback |9|, and 
a gas discharge plasma system ^^. Here we first report 
about on-off intermittency in a spatially extended dis- 
sipative system, viz. electroconvection (EC) in nematic 
liquid crystals driven by a stochastic voltage. 

EC in planarly aligned nematics is a standard system 
for pattern formation, for recent reviews see, e.g. pT| . In 
the presence of an electric field E a spontaneous fluctua- 
tion of the director leads due to the anisotropic conduc- 
tivity to a formation of space charges which tend to desta- 
bilize the homogeneously ordered state. With increasing 
strength of the driving field one observes a hierarchy of 
convection patterns of increasing complexity. The pat- 
terns depend on external parameters such as amplitude, 
frequency and wave form of the driving voltage which are 
conveniently adjustable in the experiment. The hydrody- 
namic fiow induces a modulation of the director field and 
thus of the effective indices of refraction which leads to 
transmission patterns easily observed with a microscope. 

In previous experiments, the superposition of a deter- 
ministic AC field with a stochastic field, E = E'^'^*'{t) + 
S"*"'^'* (t) , was studied. A variety of noise induced phe- 
nomena including stabilization or destabilization of the 



homogeneous state, and a change from continuous to dis- 
continuous behaviour of the threshold as a function of the 
noise strength was observed [O-nw and has stimulated 
theoretical work |l7|Jl^ ]. 

As long as the characteristic time of the noise Tstoch is 
small compared with characteristic times of the system 
the threshold towards pattern formation appears sharp 
as for deterministic driving. It is however typical for EC 
in nematics that one of the systems characteristic times 
decreases both with increasing strength of the threshold 
voltage and with increasing wave number of the pattern 
and may reach the order of Tgtoch |18| . In this case, one 
observes intermittent bursts of a regular spatial stripe 
pattern which makes a nai've experimental determination 
of the stochastic threshold difficult |l^ . A similar phe- 
nomenon was noticed in a highly doped nematic material 
for very high strength of the noise |13] where a direct 
transition towards chaos occurs via intermittent bursts 
of spatially incoherent structures embedded in a homo- 
geneous background. 

In this paper, we consider the simplest case of pure 
stochastic excitation, E — E'^*"'^'^ (t) . To achieve a sta- 
tistical characterization we have determined experimen- 
tally the distribution of the duration r of laminar, i.e. 
undistorted phases (off-states) which are interrupted by 
bursts of a stripe pattern (on-states). Approaching a crit- 
ical voltage from below, this distribution is governed by 
a power law r~^' ^ over several orders of r as it is typical 
for on-off-intermittency. This result is confirmed by sim- 
ulations of the linearized electrohydrodynamic equations 
at the sample stability threshold ||l8|] . 

We use the standard experimental set-up with a com- 
mercial cell (Linkam) providing planar anchoring of the 
director by antiparallely rubbed polyimide coatings, two 
transparent ITO electrodes of 5 x 5mm^ and a cell gap of 
d — 50/im. The nematic material is a mixture of four 
disubstituted phenylbenzoates [ |l6| which has a nematic 
range from below room temperature to 70.5°C. The cell 
temperature is controlled at 30° C by a Linkam heating 



stage. Images are recorded using a Jenapol-d polarizing 
microscope and a Hamamatsu B/W camera with con- 
troller C2400. The transmission images of light polarized 
along the director easy axis no are recorded digitally. We 
resolve the images in 500x400 (330x250) pixels with 8 
bits greyscale and calculate standard deviations in real 
time at a rate of 1/7 (1/20) s in the conductive (dielec- 
tric) regime. 

As driving process we use the dichotomous Markov 
process (DMP) which is easily generated and faciliates 
the theoretical analysis. The DMP iJ^MP j^^^pg ran- 
domly between ±E with average rate a. The distri- 
bution of times f between two consecutive jumps is 
aexp (— af); the autocorrelation decays exponentially, < 
^DMP^DMP >^ E^exp[-2a{t - t% i.e. r.toch = l/2a. 
We will refer to i' = a/2 as the mean frequency and to 
E as the amplitude of the noise. Both in experiment and 
simulation, sequences of the DMP are generated by the 
same algorithm; technically f is limited to vary between 
lO-^s and lO^s. 

For excitation by a periodic square wave, one finds a 
typical frequency dependent threshold voltage Uc {U = 
Ed) for the stability against formation of normal rolls, 
cf. Fig. |l|. There is a sharp transition at i^c = 38 Hz be- 
tween the conductive regime (oscillating space charges) 
characterized by a wavenumber fc^; « n/d and the dielec- 
tric regime (oscillating director deflections), where k^ is 
an order of magnitude larger, cf. Fig. ^ 
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FIG. 1. Threshold voltages and wave numbers for driv- 
ing with periodic and stochastic square waves of (mean) fre- 
quency V. Experimental data (diamonds: periodic case; cir- 
cles, triangles, and full squares: stochastic case, see text; 
dashed lines guide the eyes) are compared with results from 
the two-dimensional theory (full line: periodic case; dotted 
line: stochastic case). 

For stochastic driving, we have observed bursts of 
stripe pattern uniform across the system already below 
the onset threshold for a periodic voltage of the same 
frequency. The stochastic voltage has a broad frequency 
spectrum which contains low frequency contributions. 
Occasional bursts of convective pattern can be expected 
at voltages above the DC threshold (i^ = in Fig. [ip). 
With increasing voltage the frequency of the bursts in- 
creases. In Fig. |lp, we have plotted the voltages which 
correspond to a ratio of 75% (circles) and 25% (trian- 



gles) of laminar phases, respectively. The full squares at 
1/ = 60 Hz and 180 Hz indicate the voltages for which the 
experimentally determined distribution of laminar peri- 
ods T is best described by a t~^/^ law (see below). The 
theoretical results obtained from the sample stability cri- 
terion explained below agree very well with the experi- 
mental data. (Both in the periodic and the stochastic 
case we have used the same material parameters.) 




FIG. 2. Snapshots of the bursts of roll patterns (a, b) 
in the conductive regime {f = 60 Hz, k^ = 1.2 x 10'^ cm"'^, 
area size 200 x 164 /^m ) and (c, d) in the dielectric regime 
{ly = 180 Hz, k^ = 7.2 x 10^cm"\ area size 134x102 /xm^) at 
different times, (a, c) show the patterns just at the threshold 
of the on-state ((Trd ~ 0.1) whereas in (b, d) the rolls are fully 
developed ((Trd = 0.8). 

The set of images in Fig. shows the stripe pat- 
tern at different times of a burst just at the thresh- 
old of the on-state and the fully developed pattern in 
both the dielectric and the conductive regime. We char- 
acterize the intensity modulation of these patterns by 
froi = (<7 — cro)/(l — Co) where a is the normalized stan- 
dard deviation from the average intensities taken over all 
pixels of the image at a given instant, and ctq is the value 
of a for zero voltage. This procedure allows a real time 
characterization of the patterns. It is justified since the 
largest Fourier coefficient dominates the intensity modu- 
lation , and both quantities have nearly equivalent traces. 
Although the relation between cTrci and the director de- 
fiections is nonlinear, the approach used here is sufficient 
if we are interested only in the frequency and duration of 
the bursts and not primarily in their amplitudes. 

In Fig. we present trajectories of a^ciit) observed 
in the experiment for different noise amplitudes but the 
same seed of the driving DMP. As laminar phase we 
define the period in which arc\{t) is below a threshold 
Clam = 0.1 where the choice of CTiam is not crucial. The 
laminar periods are interrupted by bursts of the convec- 
tion structures, the frequency of which increases with the 
applied voltage amplitude U. 

In Figs. and a we compare the experimentally de- 
termined distributions p(t) for the occurence of laminar 
phases of duration t with those obtained in simulations 
in the conductive and the dielectric regime, respectively. 
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FIG. 3. Bursts of the intensity modulation crrei(t) in the 
conductive regime (a) just at the threshold and (b) 0.3 V 
above threshold for identical trajectories of the DMP which 
makes about 9 x 10* jumps in the period shown {ly — 60 Hz). 
The dashed lines indicate aiam above which the system is in 
the on-state and else in the off-state. 
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FIG. 4. Normalized distribution p(r) in the conductive 
regime {u — 60 Hz) just at the stability threshold and slightly 
above. Shown are results (a) from experiment for U = 14.6 
V (squares) and 15.3 V (triangles), and (b) from simulations 
for U = Uc = 18.2 V (fuU line) and U = 19.0 V (dashed line). 
Mode selection gives a wave number kx — 1484 cm~ used in 
the simulation. The dash-dotted lines indicate a power law 
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FIG. 5. Normalized distribution p(r) in the dielectric 
regime {v — 180 Hz) slightly below Uc, just at Uc, and slightly 
above. Shown are (a) experimental results for U — 76.0 V 
(circles), 79.5 V (squares), and 83.0 V (triangles), and (b) 
simulations for U = 82.0 V (dotted hne), U = Uc ^ 84.9 
V (full line), and U = 90.0 V (dashed line). Mode selection 
gives kx = 7670 cm~^ used in the simulation. Dash-dotted 
line as in Fig. 0. 



The experimental histograms contain data from about 
2500 bursts of crrci(i) recorded over 1...6 h. With increas- 
ing voltage amplitude the frequency of bursts increases, 
i.e. longer laminar periods occur less frequently. At a 
critical voltage Uc (full squares in Fig. ||a) the distribu- 
tion is governed by a power law p ~ r^/^ over several 
decades. Deviations occur for very small r due to the 
finite time resolution and for very large r due to the 
background fluctuations from thermal noise. Increasing 
the amplitude beyond the critical voltage leads to ex- 
ponential corrections to the power law discussed below. 
Increasing the voltage further, the system does not reach 
a state of permanent convection, instead the patterns be- 
come more and more spatially irregular while keeping its 
intermittent character. 

The theoretical description of the laminar phases 
is based on the linearized nemato-electrohydrodynamic 
equations. In a two-dimensional idealization using stress- 
free boundary conditions, inserting as a test mode a pe- 
riodic roll pattern characterized by wave numbers kx and 
kz = Ti/d they reduce to the ordinary differential equa- 
tions @ 



C{t)z, z- 



^ 



(1) 



where q and ip are the amplitudes of the space charge 
density and of d^^ (</? is the angle between director and 
electrode plates), and 



C(t) 



1/r, 



aE, 



DMP 



<JHEf^P 

Ai - A2E^ 



(2) 



The parameters Tg, an, a, Ai, and A2 depend on ma- 
terial properties and on the wave number k^ which is 
determined by minimizing the threshold voltage, cf. [[l8[ . 
Between two consecutive jumps of iJ^MP ^^. ^^ ^^^ t^+i 
the matrix C(i) is constant, and the time evolution is 
given by z{t) = T^" {t — t^)z{t^) for t^ < t < t^+i where 
T''(t) is the time evolution matrix for sguE^^^ — s. 
Iteration gives the formal solution [Q for a given real- 
ization of the driving process with jumps at the random 
times t^, z/ = 1, ..., n 



z{t) - T^"(i - i„) • ■ • T'«{h - to)zito) 



(3) 



The threshold voltage Uc for a given wave number is 
determined by the zero of the largest Lyapunov expo- 
nent Ai of the product of random matrices in (0) which 
can be calculated analytically as well as the second Lya- 
punov exponent A2 < Ai at the threshold |I8]. Whereas 
Tsys,i = l/|Ai| diverges at the threshold we found for 
the examples presented here in the conductive (dielec- 
tric) regime Tsys,2 = I/IA2I = 3.0 x IQ-^ s (0.9 x IQ-^ s) 
which is of the same order as Tstoch — l/2a = 4.2 x 10~'^ 
s (1.4 X 10-3 s). 

The numerical simulation generates trajectories z(t) 
starting from a small nonzero initial value z(to), cf. Eq. 



(^ . To model the background of thermal fluctuations of 
ip we introduced a lower cutoff V'min, i-e. ip ~^ V'minSgnV' 
for lipl < ■0min- In the dielectric regime we addition- 
ally reset g in a similar way. A trajectory is considered 
laminar as long as \ip\ is smaller than a given threshold 
V'lam = 2 X lO'^Vmin- At Uc the distribution is a power 
law T~^''^ over several orders of r with deviations for very 
small and very large t as in the experiment, cf. Figs. 0b 
and 0b. The range of validity of the power law increases 
when V'min/V'iam is lowcrcd. Also for voltages smaller or 
larger than Uc the shape of the simulated distributions is 
very similar to that obtained in experiment. The shoul- 
der for large r becomes more pronounced for voltages 
below the critical voltage. 
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FIG. 6. (a) Normalized distribution p(r) in the conduc- 
tive regime {u = 60 Hz) obtained from simulations with very 
small lower cutoff at Uc (full squares), and slightly above: 
AU = 0.05 V (squares) and AU = 0.2 V (triangles). The 
lines indicate the modified power law mentioned in the text, 
(b) shows the mean value of duration of laminar periods r as 
function of AU. kx as in Fig. M. 

To distinguish between effects due to thermal fluctu- 
ations and due to a deviation from the stability thresh- 
old AC/ — U ~ Uc we have also performed simulations 
with V'min = 10~^°°^iam (which is obviously too small 
for comparison with experiment). At Uc the power law 
holds now over 8 decades. Above threshold the results 
agree very well with p ^ t~^/^ exp(— const AC/^r), see 
Fig. ya. The mean duration of laminar periods behaves 
like f ^ A[/~^, cf. Fig. gb. Similar laws have been ob- 
tained analytically in g for a one-dimensional mapping. 

We have found on-off intermittency in a spatially ex- 
tended dissipative system driven by multiplicative noise 
at parameter values where the first instability is towards 
spatially regular structures. If the strength of the noise 
approaches a critical value both experiment and simu- 
lations of the electrohydrodynamic equations lead to a 
power law with exponent —3/2 for the distribution of 
laminar periods. The simulations show that this critical 
value is just the threshold of stability according to the 
sample stability criterion |18[ . 
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